Antiferromagnetic integer-spin chains in a staggered magnetic field: approaching the 
thermodynamic hmit through the infinite-size DMRG. 
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We investigate the behavior of antiferromagnetic integer-spin chains in a staggered magnetic field, by 
means of the density-matrix renormahzation group, carefully addressing the role of finite-size effects 
within the Haldane phase at small fields. In the case of spin S — 2, we determine the dependence 
of the groundstate energy and magnetization on the external field, in the thermodynamic limit, and 
show how the peculiar finite-size behavior can be connected with the crossover in the groundstate 
from a spin liquid to a polarized Neel state. 



75.10.Jm, 75.30.Cr, 75.40.Mg 
I. INTRODUCTION 

The behavior of antiferromagnetic (AFM) Heisenberg 
spin chains is dominated by quantum fluctuations which 
suppress magnetic long-range order. However, after Hal- 
dane's proposal ||l|, it has become clear that the value of 
the spin S plays a crucial role. Integer and half-integer 
spin cases are indeed completely different. The spectrum 
of the half-integer spin chain is gapless in the thermody- 
namic limit, and the system is critical, in the sense that 
the linear response to an infinitesimal staggered magnetic 
field, coupled to the would-be order parameter, is diver- 
gent (at zero temperature). Conversely, the integer-spin 
chain has a spectrum which stays gapped in the thermo- 
dynamic limit. The first excited state has a finite dis- 
tance from the groundstate (the so-called Haldane gap), 
and the system is not critical, i.e., the linear response 
to an infinitesimal staggered magnetic field is finite. The 
groundstate, which is characterized by a finite correlation 
length, is a spin liquid. 

The two behaviors are only reconciled as S" — > oo. This 
limiting case, which corresponds to the suppression of 
quantum fluctuations, leads to the closure of the Hal- 
dane gap and to the divergence of the staggered suscep- 
tibility even in the integer-spin case, which becomes thus 
indistinguishable from the half-integer-spin case, as it is 
naively expected at large S. 

All theoretical approaches are well controlled only in 
the limit of large S, while the increasing relevance of 
quantum fluctuations makes the predictions less accurate 
as S is reduced. On the other hand, the introduction of 
the staggered magnetic field gradually freezes quantum 
fluctuations, and the problem arises of the description of 
the spin-liquid state in a staggered magnetic field and 
its evolution to the frozen state. A non-linear tr-model 
approach has been recently developed to investigate this 
evolution |g] . This approach relies on the Haldane ansatz, 
and becomes exact in the limit 5 — *■ oo, in the absence 
of staggered magnetic field. However, the validity of the 
ansatz must be limited by some finite field, at which the 



spin-liquid description breaks down. Moreover, at small 
5, the starting ansatz is (quantitatively) not appropriate 
even at zero field, and corrections must be considered, 
when describing the field-driven crossover from the spin 
liquid to the polarized AFM state. 

This paper is therefore devoted to the numerical anal- 
ysis of the response of the AFM integer-spin chains to a 
staggered magnetic field, by means of the density-matrix 
renormahzation group (DMRG) |^ , which gives basically 
exact results for one-dimensional systems, that are not bi- 
ased by any a priori assumption, and represent an almost 
ideal numerical test for analytical results. 

The case S = 1 has been most extensively studied 
within DMRG [Q , but the case 5 = 2 ||] represents quite 
a hard numerical task even in the absence of a magnetic 
field, and lacks any analysis in the presence of the stag- 
gered field, that will be the main object of our numerical 
investigation. 

Our starting Hamiltonian reads 



L-l 



7^ = J^S,.S,+i-i/^(-l)*5f, 



(1) 



where S^ is the spin-S* operator on site i, S = 1,2,... 
is integer, and L is the number of sites in the chain. 
J > is the AFM Heisenberg coupling, H is the am- 
plitude of the staggered magnetic field along the z axis, 
which is coupled with the staggered magnetization M^ ~ 
^^(— l)*^^ . In view of the forthcoming DMRG study, we 
explicitly assumed open boundary conditions (OBC) in 
writing the first term in (1^). 

The Hamiltonian (1^) is manifestly not invariant un- 
der reflection with respect to the mid point of the chain 
when L is even. However, in the standard implemen- 
tation of the DMRG algorithm O], L is even, and the 
reflection symmetry is used to reduce the numerical ef- 
fort by identifying the right block with the reflection of 
the left one (for more details, see the discussion in Sec. 
P). Although there is no problem in implementing a 
non-symmetric algorithm PJTl, the Hamiltonian (nl) can 
be easily made invariant under reflection by means of 



a simple gauge transformation (i.e. a local rotation of 
the reference frame), given by (— 1)'S'? -^ Sf , and, e.g., 
{-lySf -^ Sf, while Sf coincides with Sf . Once the 
transformation is performed, the Hamiltonian ffl) is re- 
cast in the form 
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n = -jj: 
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iJ ,: iJ , 
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1=1 
(2) 



where S^ = Sf ± iSf , and the external field is now cou- 
pled to the uniform magnetization S^ = J2i^i- The 
token we pay to recover reflection symmetry is the ap- 
pearance of terms of the form S^ S^_^_i and S~S~^-f^, in 
which the spin is simultaneously raised or lowered on 
neighboring sites. 

The original Hamiltonian (nf) commutes with the z- 
component of the total spin S^ = J2i ^i ^^'^^ therefore, 
the gauge-transformed Hamiltonian (0) commutes with 
the staggered magnetization M^ = J2i{~'^y ^f ■ It is eas- 
ily realized that the groundstate of the model lies within 
the M^ = subspace at large values of the field H. At 
zero field, the groundstate is instead degenerate (in the 
presence of OBC), but, since the total spin is obviously 
integer, a state with M^ = is certainly present in the 
groundstate multiplet. We explicitly checked that (at 
least a component of) the groundstate always lies in the 
M^ = subsector for every value of the amplitude of the 
external field H. Henceforth, for the sake of definiteness, 
we assume H > 0, the case of negative H being trivially 
recovered by changing S^ -^ —S^ and M^ -^ —M^ in 
the following discussion. In particular the groundstate 
energy and magnetization are even and odd functions of 
the external field respectively. In the rest of the paper we 
refer to the properties of the model described by (||), and 
in particular to the uniform magnetization along the ex- 
ternal field, (S^). The correspondence with the original 
model (ll]) is straightforward. 

The plan of the paper is the following: in Sec. O we 
discuss the general problem of extracting the thermody- 
namic properties within the DMRG approach, in connec- 
tion with the interplay of finite-size effects and truncation 
of the Hilbert space, and we specialize the discussion to 
the case of the Haldane spin-liquid phase at small mag- 
netic field; in Sec. Ill we present our numerical results. 



mostly on the 5 = 2 case, and discuss the vario us reg imes 
which characterize the model at small (Sec. [II A ) and 



large (Sec. pIB[) magnetic field, showing how these can 
be also characterized by the peculiar form of the finite- 
size corrections to the magnetization; concluding remarks 
are found in Sec. [rv| . 

We point out that throughout the paper, whenever we 
discuss finite-size corrections, we specifically refer to cor- 
rections in the presence of OBC, which are those affecting 
our DMRG results. 



II. EXTRACTING THE THERMODYNAMIC 
PROPERTIES 

After the introduction, by S. R. White §, of the 
DMRG, a wide majority of the DMRG studies focused 
on the so-called finite-size algorithm, in which the system 
size is first increased until a chosen size is reached, and 
then stopped. Subsequently, a few further "sweeps" are 
performed to improve the basis for that system size. On 
the other hand, in the infinite-size algorithm the size is 
always increased at each step. It is clear that, for a given 
lattice size, the finite-size algorithm gives the best esti- 
mate (at fixed Hilbert space) . Conversely, if one is really 
interested in the thermodynamic limit, the information 
obtained with the infinite-size algorithm are of extreme 
value HJ^, and become crucial in the integer-spin case 
at issue in this paper, as we show below. We point out 
that extreme care has to be taken in checking the conver- 
gence of the procedure with respect to the truncation of 
the Hilbert space, before any sensible statement on the 
physical meaning of the numerical results is made. 

The algorithm we used is the standard implementation 
of the infinite-size DMRG, which is described in |g] (see 
also |3,|7||). Here we limit ourselves to a quick summary 
which serves as an introduction to the notations which 
we use in the paper. As it is customary we start with a 
4-site chain which is divided into a left and a right part 
(blocks). The density matrix of, say, the right block in 
the groundstate of the entire system (superblock) is diag- 
onalized and Nk eigenstates, corresponding to the largest 
eigenvalues, are kept as the most representative of the 
block, while the other eigenstates are truncated away, 
thus reducing the size of the Hilbert space. In the pres- 
ence of refiection symmetry with respect to the mid point 
of the chain, the left block is simply the reflection of the 
right block. Once the truncated left and right blocks are 
obtained, two more sites are added in the middle of the 
chain, and the new left and right blocks are defined, each 
as the "old truncated block + newly added site" system. 
The procedure is then iterated. After Nug iterations the 
chain contains L = 2Nb.g + 2 sites. At each step, a mea- 
sure of the reliability of the procedure is provided by the 
truncation error TZ = l~X]j=i '^i^ where Wi are the eigen- 
values of the density matrix, sorted in decreasing order. 
It follows from the definition that < 7?. < 1, and that, 
in the most general case, TZ = 1 only when all the basis 
states are kept at each iteration, unless, under special 
circumstances, a simplification occurs in the spectrum of 
the density matrix (see below. Sec. [II B). The success 



of DMRG when dealing with one-dimensional systems is 
due to their peculiar topology in the presence of short- 
range interactions. In most cases, indeed, the truncation 
error can be made small by keeping a number of states 
Nk which is much smaller than the full size of the block 
Hilbert space [(25 + 1)^/^ in the present case], and can 



be considered, for all practical purposes, independent of 
the system size L at fixed TZ. We point out, however, 
that a small value of TZ does not by itself imply a good 
convergence with respect to the truncation of the Hilbert 
space, first of all because the convergence must be di- 
rectly checked on the observed physical quantity, which 
may have a peculiar dependence on 7?,, and also because 
different observables have different rates of convergence 

@- 

The observables are extracted by means of the stan- 
dard DMRG procedure H . In this paper we mainly dis- 
cuss the groundstate energy 

Eo{J,H-L,Nu) = (\n\' 

and magnetization 

M{J, H; L, Nk) = / [S' 

where the symbol (•) denotes the expectation value over 
the DMRG approximate groundstate wavefunction, for 
the system of size L, with Nk states retained in the 
density-matrix truncation procedure, and the superscript 
T means that the corresponding operator has been prop- 
erly truncated onto the superblock basis. 

In the following we define the groundstate energy per 
site in units of J, eo = Eq/{JL), and the magnetization 
per site in units of S*, m = M/{SL), such that < to < 1. 
We also introduce the notation h = H/J, since eo and m 
depend on J and H only through their ratio h. 

A glance at the magnetization curves as a function of 
the system size L, for various magnetic fields, and a fixed 
number of states (e.g., Nk = 50, see Fig. pi), shows that 
there are two different behaviors at small and large field 
respectively. The convergence of the magnetization to 
the thermodynamic limit is non-monotonic at small h, 
and m is characterized by a well-pronounced bump at a 
characteristic length Lb which increases with increasing S 
and decreases with increasing h. This property strongly 
suggests that system sizes smaller than L^ are definitely 
not representative of the thermodynamic limit, and that 
only for L ^ Lb, where a monotonic convergence to the 
thermodynamic limit is recovered, the finite system is 
a reasonable representation of the infinite-length chain. 
For 5 = 1, at the smallest field we investigated, h = 
0.002, we found Lb ~ 25 (see lowest curve in the left panel 
of Fig. p. This effect is more dramatic as S increases 
from 1 to 2. Fig. H^ clearly indicates that, for 5 = 2, 
and for the smallest applied field h = 0.0001, Lb ~ 70, 
and that, even for L ~ 200, the effect of the bump is 
not negligible, in the sense that the system appears to be 
far from an asymptotic convergence (see lowest curve in 
the right panel of Fig. m . We shall address the physics 



underlying this behavior in Sec. III. 

It should be clear that the presence of the bump 
strongly limits the possibility to draw reliable conclusions 



from finite-size DMRG calculations, unless the latter are 
pushed to exceedingly large sizes, with huge numerical 
effort. On the other hand, only a careful extrapolation 
procedure can give insights on the actual thermodynamic 
limit. The infinite-size algorithm is certainly suitable for 
this kind of study. 

Within our procedure we extract the thermodynamic 
limit for any given number of states Nk, by fitting the 
large-L behavior of the energy and magnetization as 
eo{h-L,Nk) = e'^{h-Nk) + A/L and m{h;L,Nk) = 
m°°{h;Nk) + B/L respectively. The long-living correc- 
tions are 0(1/L) due to the OBC. We explicitly checked 
that the linear fit in 1/L at small field is accurate only 
if L :^ Lb- For the smallest field (and S = 2), we had 
to consider L ~ 1600 to obtain a reliable result, and the 
situation gets worse at larger S Q. 

Once eg°(/i; Nk) and m°°{h; Nk) are obtained, we care- 
fully check the convergence with respect to the Hilbert 
space, for both m and eo, by taking the extrapolation 
to the limit Nk — > c» in the region where the observ- 
ables display a linear dependence on 1/Nk |§|,0. Up to 
Nk = 150 states were needed in the case S* = 2 for the 
smallest field h = 0.0001, to obtain the linear scaling. 
We emphasize that this procedure explicitly deals with 
the convergence to the limit of infinite Hilbert space of 
the expectation values of the different observables, and is 
not equivalent to fixing a small truncation error a priori. 



III. RESULTS AND DISCUSSION 

In this section we mostly devote our analysis to the 
case 5 = 2, for which no results in the presence of the 
magnetic field can be found in the literature. 

The groundstate energy per site eo extrapolated to 
the infinite-size limit (L — > cx)) and to the full Hilbert 
space {Nk — > cxd), is a smooth monotonically decreasing 
function of the magnetic field, correctly interpolating be- 
tween the zero-field value eo(/i = 0) = —4.76126(1), and 
the large-field asymptotic behavior eo(ft. ^ 1) ~ —Sh. 
The small-field region is plotted in Fig. g, and the 
large-field region is shown in Fig. y. Our result for 
eo(/i = 0) is obtained by an extrapolation as a function of 
h, and the (relatively) large error bar is due to this pro- 
cedure. The value compares reasonably well with Ref. 
lH (eo = -4.761244(1)). Our value is sfightly lower, sug- 
gesting that our extrapolation procedure is slightly more 
efficent that the one in Ref. B. 

The magnetization along the magnetic field can be ei- 
ther computed as outlined in Sec. O, by directly evaluat- 
ing the expectation value of the corresponding operator 
S^ on the (approximate) DMRG groundstate wavefunc- 
tion, or exploiting the Hellmann-Feynman theorem. The 
latter implies {Sz) = —dEo/dH when Eq is an exact 
eigenvalue, m may therefore be estimated by a numeri- 
cal derivative of the groundstate energy per site eo with 



respect to the field ampHtude h. It must be pointed out 
that when the wavefunction is not an exact eigenstate, 
as in the case of DMRG, the equahty does not (necessar- 
ily) hold. However, if the state obtained within DMRG 
is close to an exact eigenstate, the Hellmann-Feynman 
relation is expected to hold within good accuracy. 

The magnetization values obtained following the two 
procedures are shown in Fig. ^ and fall on a single reg- 
ular curve. This remarkable agreement represents then 
a strong evidence of the convergence of our method with 
respect to the truncation of the Hilbert space, that may 
be used as a further check, besides the control of the 
truncation error TZ, which is standard in DMRG p] , and 
the extrapolation to N^ — > oo explained in Sec. ^. 

The magnetization is a smooth function of h for all 
the values of the magnetic field, and correctly repro- 
duces both the zero-field limit m = 0, and the /i ^ 1 
limit (shown in Fig. g), where the curve saturates to 1. 
Nonetheless, a relatively sharp crossover at /i ~ 0.005 
is signaled in Fig. |j by a rapid change in the slope 
of the curve. The small-h region is characterized by a 
spin-liquid behavior, and the crossover can be naturally 
associated with a change in the nature of the ground- 
state from the spin-liquid state to a polarized Neel state, 
induced by a gradual freezing of quantum fluctuations. 
A further support of this interpretation comes from the 
analysis of the curves of to as a function of L shown in 
Fig. n. The small-ft, region coincides indeed with the 
range in which the bump is observed in the curves. This 
finding sheds light on the physical origin of the bump. 
The bump can in fact be associated with a spin-liquid 
behavior, since it is observed only in this region. 

If we now turn to the magnetization curves of Fig. n, 
we can interpret the non-monotonicity in a simple way. 
When the system has a spin-liquid groundstate, it is char- 
acterized by a spin gap, and the correlation functions de- 
cay with a characteristic lengthscale ^. If L <C £,, the 
system is too small to display the spin-liquid features, 
and behaves as if ^ ^ oo, i.e., similarly to a half-integer 
spin system, which responds more strongly to a small 
staggered external field. As a result, this leads to an 
"overshooting" of the magnetization, that assumes val- 
ues larger than the thermodynamic limit. Only when 
L ^ ^, the finite-system is representative of a spin liq- 
uid: the screening of the field becomes effective and the 
magnetization relaxes to its bulk value. This leads to an 
interpretation of Lb as an intrinsic lengthscale related to 
some correlation length ^. A more quantitative study of 
this relationship requires a model for the evolution of the 
Haldane spin-liquid phase in the small-field region, and 
is currently underway [QOl. 



A. The small-field region 

As mentioned above, it is well known that the integer- 
spin AFM Heisenberg chain has a spin-liquid ground- 
state, characterized by a spin gap, the so-called Haldane 
gap. It is quite natural that the spin-liquid behavior ex- 
tends to small finite fields. As a consequence, at the 
smallest field h ~ 0.0001, the magnetization is linear in 
the field, as opposed to the half-integer spin case, which is 
critical (i.e. with divergent linear susceptibility) at zero 
field. Our (lower-bound) estimate of the linear suscepti- 
bility is x{h = 0) = ^\h=o = 1070, in extremely good 
agreement with the calculations in Ref. |2]. In this re- 
gion the non-linear cr-model is therefore the most correct 
field theory. Unfortunately, the Haldane gap decreases by 
increasing the spin S, so that the field-driven crossover 
from the spin liquid to the antiferromagnet is expected 
to occur at an external field which decreases with in- 
creasing S (he - 0.08 for S" = 1 and h^ ~ 0.005, for 
S" = 2, as stated above). For h > he, the discrepancy 
with the non-linear cr-model results increases, even if a 
broad intermediate region exists, where the system is no 
longer a spin liquid, but the groundstate is not yet well 
described as a perturbation of the fully polarized antifer- 
romagnet. In this intermediate region, the convergence of 
the magnetization to the thermodynamic limit is mono- 
tonic, and the asymptotic value is approached from be- 
low, contrary to the spin-liquid regime. Finally, a much 
smoother crossover takes place at larger fields, signaling 
the onset of a large-field "perturbative" region, where m 
monotonically approaches its limiting value from above, 
as we discuss in Sec. 
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We point out that the field- 
driven freezing of the quantum fluctuations is seen in the 
gradual simplification of the spectrum of the density ma- 
trix, i.e., a continuous reduction of the truncation error 
TZ by increasing h at fixed Nk- 



B. The large- field region 

In the limit h —^ oo the groundstate is trivially polar- 
ized. In this limit the spectrum of the density matrix 
is therefore dominated by a single state, with eigenvalue 
wi = 1, and the truncation error TZ vanishes as soon 
as Nk > 1. It is therefore reasonable to expect that 
at large magnetic fields, h ^ 1, the spectrum of the 
density matrix will be much simpler than in the spin- 
liquid phase, making the convergence to the full Hilbert 
space relatively easy. As a consequence, the DMRG 
data are exceedingly accurate a.t h = 1, even by keep- 
ing Nk ~ 50 states (giving 7^ -- 10""), and at h = 10, 
keeping Nk — 30 states (giving TZ ^ 10~^^). Moreover, 
in this limit we can check our numerical data by compar- 
ison with the analytical results of standard perturbation 
theory around the J = fully polarized groundstate. 



The groundstate wavevector for the Hamiltonian (g) at 
J = is |0) = \S)i...\S)l, with a groundstate energy 



S, 



Q -^ = —HSL. As the perturbation Hamiltonian is 

L-l 



the groundstate energy up to second order in J is given 

by 



E{) ~ E. 



(0) 



Vn 



E' 



iVfeoP 



p(0) p(0) ■ 



The excited states which contribute to the sum over k are 
all degenerate, and correspond to wavevectors of the form 
|fc) ^\S)i...\S^ l)k\S-l)k+i ■ ■ . \S)l, k^l, ...,L - 1, 
with unperturbed eigenvalues Ef. = —HSL + 2H. As a 
result, for /i ^ 1 we find 



eo 



'hS 



1 + S 



S^ 



1 
2h 



(3) 



where the infinite-size value (first term in the r.h.s.) and 
the leading finite-size corrections 0(1/L) (second term in 
the r.h.s.), are explicitly shown. The agreement with the 
DMRG data for ft, > 1 is perfect, as shown in Fig. 0. 

Analogously one can obtain the magnetization along 
the field, by calculating the groundstate wavefunction up 
to second order in perturbation theory, and then evalu- 
ating the expectation value of the operator S^ . It is how- 
ever much simpler to make use of the Hellmann-Feynman 
theorem once more and calculate {S^) = ~dE^/dH ~ 
SL - S^{L - l)/(2/i2), obtaining 



S 
2ft2 



S 



2h^L' 



(4) 



We have thus found that the convergence to the thermo- 
dynamic limit [the first term in the r.h.s. of Eq. (j^] is 
monotonic and the finite-size corrections [the second term 
in the r.h.s. of Eq. (^) are positive, in perfect agreement 
with the DMRG data (see Fig. |^. The deviation from 
the saturation value in the thermodynamic limit vanishes 
as l//i^, in disagreement with the l/v" behavior found 
within the nonlinear a- model approach ^. Since the 
calculation of Ref. g] relies on the Haldanc ansatz, the 
disagreement is expected at large magnetic fields, where 
the groundstate is Neel-like, rather than spin-liquid-like. 



the finite-size version, to extract the thermodynamic 
properties of these systems, due to the presence of size- 
able characteristic lengthscales, which induce anomalous 
finite-size effects. 

We have been able to determine three regions in the 
phase diagram of the S — 2 chain, as a function of the ap- 
plied staggered field. The first region, for < /i ^ 0.005 
is the spin-liquid region, in which the spectrum is char- 
acterized by the Haldane gap, the spin susceptibility is 
linear in the field, with a coefficient x(0) — 1070, that 
compares extremely well with a non-linear cr-model cal- 
culation 0] . In this region the convergence of the mag- 
netization to the thermodynamic limit is not monotonic 
as the system size L is increased, due to the existence of 
a characteristic lengthscale ^. Thus, quite large systems 
must be considered to correctly describe the spin-liquid 
phase. Indeed, for L <C f the system is expected to be- 
have as if ^ — !■ CX3, as it is the case for half- integer spin, 
with a strong increase of the magnetization. It is only for 
i ^ ^ that the screening of the external field in the spin- 
liquid phase becomes effective and the magnetization can 
relax to its bulk value. 

The second region (0.005 ^ h < 1) is an intermediate 
crossover region, where the system is no longer a spin 
liquid, due to the field-driven freezing of quantum fluctu- 
ations, but the field is not large enough to give rise to a 
linear antiferromagnet. As a result, both the non-linear 
cr-model (valid in the Haldane phase), and perturbation 
theory around the linear antiferromagnet (obviously valid 
for h —^ cxj), are not accurate. In this region the magne- 
tization curve monotonically approaches the thermody- 
namic limit, and the limit is reached from below (contrary 
to the Haldane phase). 

Finally, for h > 1, the magnetic field is large enough 
to completely freeze quantum fiuctuations and give rise 
to an almost saturated Neel antiferromagnet, that can 
be well described by perturbation theory in l/h. As 
predicted by perturbation theory, the magnetization de- 
creases monotonically with increasing system size and ap- 
proaches its thermodynamic limit from above. 

A more detailed and quantitative comparison between 
our numerical results and the field-theoretical results of 
the nonlinear cr-model is presently in progress, aiming 
to clarify the identification of the characteristic length 
Lfc observed in the DMRG data with some correlation 
length, and the precise region of validity of the analytical 
approaches [M. 



IV. CONCLUSIONS 

We have presented extensive DMRG calculations for 
integer-spin antiferromagnetic Heisenberg chains in an 
external staggered magnetic field. The infinite-size 
DMRG algorithm has proven much more suitable than 
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FIG. 1. Average magnetization jn as a function of the 
length of the chain L, in the case S = 1 (left) and S = 2 
(right), for Nk = 50. Notice the different scale on the L-axis. 
The different curves correspond to increasing values of the 
magnetic field, from bottom to top: in the left panel h = 
0.002, 0.004, 0.006, 0.008, 0.01, 0.02, 0.03, 0.04, 0.06, 0.08, 
0.1, 0.12, 0.14, 0.16, 0.18, 0.2, 0.3, 0.4, 0.5, 0.75, 1.0, 2.0; 
in the right panel h = 0.0001, 0.0003, 0.0006, 0.0009, 0.002, 
0.004, 0.006, 0.008, 0.01, 0.02, 0.03, 0.04, 0.05, 0.1, 0.2, 0.5, 
1.0, 2.0. 
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FIG. 2. Groundstate energy per site in the case S = 2, 
extrapolated to the thermodynamic limit L — > oo and to 
Nk -^ cxD (full circles), as a function of the magnetic field. 
The solid line is a guide to the eye. 
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FIG. 3. Large-field behavior of the groundstate energy per 
site in the case 5* = 2 (full circles) compared with the sec- 
ond-order perturbation theory result, Eq. (h) for L — > oo 
(dashed line). To make the comparison more transparent 
~eo{h)/Sh — 1 is plotted against 1/h, in log-log scale. 
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FIG. 4. Average magnetization per site m in the case 
S = 2, extrapolated to the thermodynamic limit i — > cx3 
and to Nk — > C)0, as a function of the magnetic field. The 
full circles are direct DMRG measures, joined by the solid 
line as a guide to the eye, while the open dots are obtained 
by numerically differentiating the energy curve (see Fig. H) 
according to the Hellmann-Feynman theorem. 
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FIG. 5. Large-field behavior of the average magnetization 
per site m in the case S = 2 (full circles) compared with the 
perturbation theory result, Eq. M) for L ^ oo (dashed line). 
To make the comparison more transparent 1 — m is plotted 
against 1/h, in log- log scale. 
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